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Abstract 

We investigate spin-dependent parton distributions in the polarized 
virtual photon target in perturbative QCD up to the next-to-leading 
order (NLO). In the case < < Q^, where -Q^ (--P^) is the 
mass squared of the probe (target) photon, parton distributions can 
be predicted completely up to NLO, but they are factorization-scheme- 
dependent. We analyze parton distributions in six different factorization 
schemes and discuss their scheme dependence. We study, in particular, 
the QCD and QED axial anomaly effects on the first moments of par- 
ton distributions to see the interplay between the axial anomalies and 
factorization schemes. We also show that the factorization-scheme de- 
pendence is characterized by the large-x behaviors of quark distributions 
in the virtual photon. Gluon distribution is predicted to be the same up 
to NLO among the six factorization schemes examined. In particular, 
the first moment of gluon distribution is found to be factorization-scheme 
independent up to NLO. 
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1 Introduction 



In the two-photon process of e~^e~ colhsion experiments, we can measure the struc- 
ture functions of the virtual photon (Fig.l). The advantage in studying the virtual 
photon target is that, in the case 

A2<p2<g2 (1.1) 

where —Q^ {~P'^) is the mass squared of the probe (target) photon, and A is the 
QCD scale parameter, we can calculate the whole structure function up to the next- 
to-leading order (NLO) in QCD by the perturbative method, in contrast to the case 
of the real photon target where in NLO there exist non-perturbative pieces [0, §]. 
The spin-independent structure functions F2'(x, Q^,P^) and F2(x, Q^, -P^) as well 
as the parton contents were studied in the leading order (LO) and in NLO 0-111 • 
The target mass effect of unpolarized and polarized virtual photon structure in LO 
was discussed in Ref . p!0[] . 

The information on the spin structure of the photon would be provided by the 
resolved photon process in polarized version of the DESY electron and proton collider 
HERA |ri], More directly, polarized photon structure function can be measured 
by the polarized e^e~ collision in the future linear colliders. For the real photon 
(P^ = 0) target, there exists only one spin- dependent struture function, gJi^XjQ"^), 
which is equivalent to the structure function W7(x, Q^) {gj = 2W2) discussed some 
time ago in [0, [1^. The LO QCD corrections to gj for the real photon target was 
first calculated by one of the authors ||15| and later in Refs.|[l^, [l^, while the NLO 



QCD analysis was performed by Stratmann and Vogelsang |18[. The first moment 
of the photon structure function gl has recently attracted attention in the literature 



[|T^, [19], ^ in connection with its relevance for the axial anomaly. More 



recently the present authors investigated ||23| the spin-dependent structure function 



gi{x, Q^, P^) of the virtual photon up to NLO in QCD, where P^ is in the above 
kinematical region ( p..l| ). The analysis was made in the framework of the operator 
product expansion (OPE) supplemented by the renormalization group method and 
also in the framework of the QCD improved parton model using the DGLAP 
parton evolution equations. 
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In the past few years, accuracy of the experimental data on the spin-dependent 
structure function gi of the nucleon has been significantly improved . Using these 
experimental data together with the already existing world data, several groups [p6| - 
[p9[] have carried out the NLO QCD analyses on the polarized parton distributions 
in the nucleon. These parton distributions may be used for predicting the behaviors 
of other processes such as polarized Drell-Yan reactions and polarized semi-inclusive 
deep inelastic scatterings, and etc. However, parton distributions obtained from the 
NLO analyses are dependent on the factorization scheme employed. It is possible 
that parton distributions obtained in one scheme may be more appropriate to use 
than those in other schemes. In the case of nucleon target, however, it may be 
difficult to examine the features of each factorization scheme, since for the moment 
it is inevitable to resort to some assumptions in order to extract parton distributions 
from the experimental data. 

On the other hand, it is remarkable that, in the case of virtual photon target 
with a virtual mass —P^ being in the kinematical region of Eg. ( |1.1| ), not only the 
photon structure functions but also the parton distributions in the target can be 
predicted entirely up to NLO in QCD. Thus, comparing the parton distributions 
predicted by one scheme with those by other schemes, we can easily examine the 
features of each factorization scheme. In consequence, the virtual photon target may 
serve as an optimal place to study the behaviors of parton distributions and their 
factorization-scheme dependence. 

In this paper we examine in detail the polarized parton (i.e., quark and gluon) dis- 
tributions in the virtual photon target. The polarized parton distributions are partic- 
ularly interesting due to the fact that they have relevance to the axial anomaly fS^ . 
The interplay between the QCD axial anomaly and factorization schemes has already 
been discussed for the spin-dependent structure function gi of the nucleon ||3T[|-||37[ 



It was explained there that the QCD axial anomaly effect is retained in the fiavor- 
singlet quark distribution in the nucleon in the standard MS scheme, but it is shifted 
to the gluon coefficient function in such a scheme called chirally-invariant (CI) fac- 
torization scheme. Now it should be pointed out that the polarized photon target 
is unique in the sense that not only QCD but also the OED axial anomaly takes 
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place. The QED axial anomaly, which is U(l) anomaly, emerges when a quark has 
an electromagnetic charge. Thus the flavor-non-singlet quark distribution is also 
relevant, besides the flavor- singlet one. Depending upon factorization schemes, the 
QED axial anomaly effect resides in both the flavor-singlet and non-singlet polarized 
quark distributions in the virtual photon, or it is shifted to the photon coefficient 
function, in which case we arrive at an interesting result: the flrst moments of the 
polarized quark distributions in the virtual photon, both flavor singlet and non- 
singlet, vanish in NLO. Also we flnd that the large a;-behaviors of polarized quark 
distributons dramatically vary from one factorization scheme to another. Indeed, 
for x — i> 1, the quark distributions positively diverge or negatively diverge or remain 
flnite, depending on factorization schemes. 

We perform our analyses in six different factorization schemes, (i) MS, (ii) CI 
(chirally invariant) (it is also called as JET) [^, |3^, (iii) AB (Adler-Bardeen) ||36[| , 
(iv) OS (off-shell) ||, (v) AR (Altarelli-Ross) [||, and finally (vi) DIS^ schemes [||], 
and see how the parton distributions change in each scheme. In particular, we study 
the axial anomaly effects on the first moments and the large-x behaviors of parton 
distributions. Gluon distribution in the virtual photon is found to be the same up 
to NLO, at least among the factorization schemes considered in this paper. Further- 
more, the first moment of gluon distribution turns out to be factorization-scheme 



independent up to NLO. Part of the result has been briefiy reported elsewhere |^ 



In the next section we discuss the polarized parton distributions in the virtual 
photon. The explicit expressions for the fiavor singlet- (non-singlet-) quark and gluon 
distributions predicted in QCD up to NLO are given in Appendix A. In Sec. 3, we 
derive the transformation rules for the relevant two-loop anomalous dimensions and 
one-loop photon matrix elements from the MS scheme to other factorization schemes 
and then explain particular factorization schemes we consider in this paper. In 
Sec. 4, we examine the first moments of parton distributions with emphasis on the 
interplay between the QCD and QED axial anomalies and the facorization schemes. 
The behaviors of parton distributions near x = 1 and their factorization-scheme 
dependence are discussed in Sec. 5. The numerical analyses of parton distributions 
predicted by different factorization schemes will be given in Sec. 6. The final section 
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is devoted to the conclusion and discussion. 

2 Polarized parton distributions in photon 

Let q].{x,Q'^,P'^), CJ-lx , Q"^ , P^) , r].{x , Q"^ , P"^) be quark with z-flavor, gluon, and 
photon distribution functions with ± hehcities of the longitudinally polarized virtual 
photon with mass —P^. Then the spin-dependent parton distributions are defined 
as Aq' = qi + q'_^-ql-ql , AG^ = G\-Gl , and AV^ = V\-Vl . In the leading 
order of the electromagnetic coupling constant, a = e'^/An, AV^ does not evolve 
with and is set to be Ar'''(a;, Q^, P^) = 5(1 — x). For later convenience we use, 
instead of Ag*, the flavor singlet and non-singlet combinations of spin-dependent 
quark distributions as follows: 

Aq-^s = E 

i 

Mis - Ee?K-^)- (2-1) 

In terms of these parton distributions, the polarized virtual photon structure func- 
tion gi{x, Q^, P^) is expressed in the QCD improved parton model as 

gj{x, Q^ P') = /' ^ {Aq^iy, Q^ P') AC](-, Q') + AG\y, Q\ P') AC^-, Q') 

y ^ y y 

+Aqlsiy,Q',P') AClsi^,Q')] + AC:;ix,Q') . (2.2) 

where ACg{AClg), ACq, and AC^ are the coefficient functions corresponding to 
singlet (non-singlet)-quark, gluon, and photon, respectively, and are independent of 
P^. The Mellin moments of gj is written as 

gUn, Q\ P') = AC\n, Q') ■ Ag^(n, Q^ P') , (2.3) 

where 

AC"(n,Q2) = (AC] , AC^a , ^C],s , Aq) , 
Aq^{n,Q',P') = [Aq^s , , Mis , AF^) , 
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and the matrix notation is implicit. 

The exphcit expressions of Ag^, AG''', and Aqjfs ^'^ NLO can be derived 
from Eq.(4.46) of Ref.[^, which are given in Appendix A. They are written [^] 



in terms of one-(two-) loop hadronic anomalous dimensions A7°'"' (A7IJ'''") = 
ip,G) and A^y'^^ {Ins"^)^ one-(two-) loop anomalous dimensions AK^'^' (Afcf^'''") 
{i = ipjG, NS) which represent the mixing between photon and three hadronic 
operators R" {i = tp, G, NS), and finally AA", the one-loop photon matrix elements 
of hadronic operators renormalized at /x^ = P^(= —p^), 

ilip) I R-if^) I 7(P))U-P^ = ^AA- (^ = ij,G,NS) . (2.4) 

The photon matrix elements AAf are scheme- dependent. In one- loop order, they 
are given, in the MS scheme, by |34 



i^, MS ^g4^ _ ^g2^2 "-^^NS, MS 



12(e^)A^/|-^^5iH + r4n^-^--[ ' (2-5) 
\ n{n + 1) (n + 1)'^ n\ 



^"^G, MS - 



where Si {n) = I]?=i ^ 



3 Factorization schemes 



3.1 Transformation rules from MS scheme to a-scheme 

Although gj is a physical quantity and thus unique, there remains a freedom in the 
factorization of gJ into AC"' and Aq^. Given the formula Eq.( p.3| ), we can always 
redefine AC'^ and Aq^ as follows P2[] : 

AC\n,Q')^AC\n,Q')l = AC\n,Q^)Z-\n,Q^) , (3.1) 
Aq^n,Q\P^)^Aqin,Q\P^)\a = Za{n,Q^) Aq^{n,Q\P'') , (3.2) 

where AC^\a and Aq\a correspond to the quantities in a new factorization scheme-a. 
Note that the coefficient functions and anomalous dimensions are closely connected 
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under factorization. We will study the factorization scheme dependence of parton 
distribution up to NLO, by which we mean that a scheme transformation for the 
coefficient functions is considered up to the one-loop order, since a NLO prediction 
for gj is given by the one-loop coefficient functions and anomalous dimensions up 
to the two-loop order. 

The most general form of a transformation for the coefficient functions in one- 
loop order, from the MS scheme to a new factorization scheme-a, is given by 

ACi: = ACJ!^ -(e^)|A^,(n,a), 

AC^a\: = ^C^cm -(e^)^A.(n,a), 

AC]^s:a - AC;,"^ -^A^^5(n,a), (3.3) 

^^^T - -^3(e^)Ai(n,a), 

where (e^) = Y.i^t/Nf. The flavor-singlet (nonsinglet) quark coefficient functions 
are expanded up to the one-loop order as 

ACl'- = {e%l + ^AB-s + 0{al)) , (3.4) 



47r 

A 

47r 

with ASg = AS]^5 . The Az{n, a) (Ai(n, a)) term tells how much of the QCD 
(QED) axial anomaly effect is transferred to the coefficient function in the new 
factorization scheme. The gluon and photon coefficient functions ACq " and AC^' " 
start from the one-loop order (i.e., from the NLO): 

ACS'" = {e'){^ABl + 0{al)) , 

Aq'" = ^37V^(e^)(AS; + 0(a,)) . (3.5) 

In the MS scheme, AC^' ^ has been obtained from AC^' with changes: o.sj'l'n 
{2a /as) X (Q!s/27r) , (e^) 3(6"^), and 3 is the number of colors. Thus we have 

MS = j-A^S, MS • (3-6) 



6 



Since, in the leading order, coefficient functions are given by 



AC^|Lo = AC]|Lo = ((e'), 0,1,0), 



(3.7) 



the relations ( p.3|) between the coefficient functions in the a-scheme and MS scheme 
lead to Z~^(n, Q^), which is expressed as 



fl^Aws{n,a) ^Azin,a) p{e^)Az{n,a) \ 



^AwNsin,a) f 3((e4) - (e2)2)Ai(n, a) 

V J 

(3 



where / is a 4 x 4 unit matrix. 

Now we derive corresponding transformation rules from MS scheme to a-scheme 
for the relevant two-loop anomalous dimensions. The parton distribution functions 
Aq'{n, Q^, -P^) satisfy the following evolution equation [||, ||, |5[: 

dAq\n,Q^,P^^ 



AP{n,Q'')Aq\n,Q\P^^ 



(3.9) 



where 



AP(n,g' 



fAR^4n,Q') AP^Gin,Q') Aksin,Q') \ 

APG^(n,Q2) APGG{n,Q') AkG{n,Q') 

APNs{n,Q^) AkNs{n,Q^ 





V 

From Eq. (|3.2|) , we obtain 
dAg^(n,Q^p2)|^ dZ,{n,Q^) 



J 

(3.10) 



.2 o2M , V ^2,dAq^in,Q\P^)\^ 



Aq''{n,Q\P')\^+Zain,Q' 



AP{n,Q')UAq'^{n,Q\P')\a 



d lnQ2 



(3.11) 



with 



AP(n,Q' 



dlnQ2 



+ Za{n,Q^)AP{n,Q^) 



MS 



Z-\n,Q'). (3.12) 
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The sphtting functions APi{n, Q^) {i = i/ji/j, ^jG, Gijj, GO, and NS) and Akj{n, Q^) 
{j = S, G, NS) are expanded as 



27r 



2n 



Since the QCD effective couphng constant ^^(Q^) satisfies 



(3.13) 

(3.14) 

(3.15) 



where /3o = 11 — |A^/ is the one-loop coefficient of the QCD beta function, and the 
n-th anomalous dimensions are defined as 



Akf\n) 
we find for one-loop 



--A7°'" , 
4 ■' 



APi^^(n) = -^A#- 



Ak^l\n) = -AK) 
J 8 



(l),n 



A7°''^ = A7°''^ 
a 'i, MS ' 



Aii-"'" = Ai^°'" 



j, MS ' 



(3.16) 
(3.17) 

(3.18) 



and for two-loop 

A7i/)i/), a 



^IGG, a 

A<)'r 



A' 



''l^i'^MS + 2A^(n, a) A7g'^ + 4/3oAws(n, a) 



A7^'^' 

'wG, MS 



Air 



(l),n 



■ 2A;^(n,a) [A7S'S - A7;'; + 2/3o; 
+2A«;5(n,a)A7jg , 

^7S"ms-2A^^(^'«)^7^';' 
^^GcrMs ~ 2A^(n, a) A7^^ , 
^I'^i^^j^ + ^PoAwNsin^a) , 

(l)£l 

5, M 

AX«'^ + 4Ai(n,a)3(e^)A7S, 
^<5rMs + 2A^^^5(n,a)Aif^5 

+4Ai(r^,a)3((e^)-(e^)2)A7)(f^ 



(3.19) 



Ai^l^'f^ + 2A^5(n, a) AX°'" + 4Az{n, a)3{e')Aji; , 
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The one-loop photon matrix elements of the hadronic operators, AA^ and AA^^ 
in Eq.(|2.4|), are related to each other as 



AA-^s = ' ' , (3.20) 



(A^' 7- + {e')AA; + AA^s) (3-21) 



and the sum 

a 

is factorization-scheme-independent in one-loop order [^. Thus we obtain from 



AAl, = AAl^ + 12{e')Az{n,a) , 

AAS,„ = AAI^^ = 0, (3.22) 
AA-^s,a = AAl^^^ + 12i{e')-{eY) Azin,a) . 

Note that AAq = in one-loop order. 

It is possible to choose Az{n, a) and Az{n, a) arbitrarily. In the following, we 
take Az{n, a) = Az{n, a) in the Cl-like schemes and Az{n, DIS^) ^ Az{n, DIS^) = 
in the DIS^ scheme. In one-loop order we have Aws{n, a) = AwNs{n, a). Thus from 
now on, we set Aws{n, a) = AwNs{n, a) = Aw{n, a). Let us now discuss the features 
of several factorization schemes. 



3.2 The MS scheme 

This is the only scheme in which both relevant one-loop coefficient functions and 
two- loop anomalous dimensions were actually calculated ||3^ , In fact 



there still remain ambuguities in the MS scheme, depending on how to handle 75 
in n dimensions. The MS scheme we call here is the one due to Mertig and van 
Neerven and Vogelsang |^8|] , in which the first moment of the non-singlet quark 



operator vanishes, corresponding to the conservation of the non-singlet axial cur- 
rent. Indeed we have A7^]:^'"^ = 0. Explicit expressions of the relevant one- loop 
coefficient functions and two-loop anomalous dimensions can be found, for example, 
in Appendix of Ref. It is noted that, in the MS scheme, both the QCD and 
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QED axial anomalies reside in the quark distributions and not in the gluon and 
photon coefficient functions. In fact we observe 

= 24C^^/^0' (3-23) 
^B-h-s = ^^rMs = 0- (3.24) 

where Cp = ^ and Tf = Also we find from Eg. ( |2.5| ) that the ffist moments of 
the one-loop photon matrix elements of quark operators gain the non-zero values, 
i.e., 

A^rk = J^^yz\^2 ^^Ts, MS = -12(e^)iV, , (3.25) 
which is due to the QED axial anomaly. 

3.3 The Cl-like schemes 



The EMC measurement of the ffist moment of the proton spin structure function 
QilXjQ"^) presented us with an issue called "proton spin crisis". Since then many 
ideas have been proposed as solutions. One simple and plausible explanation was 
that there exists an anomalous gluon contribution to the ffist moment [pT|-p3[ origi- 



nating from the QCD axial anomaly. This explanation was later supported with 
a notion of the factorization-scheme dependence. There is a set of the factorization 
schemes in which we obtain 

AB^r' = -2iV;, A7S'"=^ = . (3.26) 

Let us call them Cl-like schemes. In this paper we consider four Cl-like schemes, in 
which we take Az{n,a) = Az{n,a), since both QCD and QED anomalies originate 
from the similar triangle diagrams. With this choice, the relation between the one- 
loop gluon and photon coefficient functions, which holds in the MS scheme, also 
holds in the Cl-like schemes, 

ci-iike = -^^^Bl ci-iike • (3.27) 



Thus, in addition to the relations in Eg. ( |3.26D , we obtain in the Cl-like schemes 



A-B; Cl-like — ~4 , AA^^ Cl-like — ^^NS, Cl-like " . (3.28) 
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(i) [The chirally invariant (CI) scheme] In this scheme the factorization of the 
photon-gluon (photon-photon) cross section into the hard and soft parts is made so 
that chiral symmetry is respected [p7| , p8[| and the QCD and QED anomaly effects are 
absorbed into the gluon and photon coefficient functions. Thus the spin- dependent 
quark distributions in the CI scheme are anomaly-free. The transformation from 
the MS scheme to the CI scheme is achieved by 

Aw{n, a = CI) = , Az{n, a = CI) = Az{n, a = CI) = ^Nf ^^^^ . (3.29) 



It has been argued by Cheng and Miiller and Teryaev that the x-dependence 
of the axial-anomaly effect is uniquely fixed and that its x-behavior leads to the 
transformation rule ( |3.29| ) and thus to the CI scheme. 

(ii) [The Adler-Bardeen (AB) scheme] Ball, Forte and Ridolfi [|S^ proposed several 
Cl-like schemes for the analysis of the nucleon spin structure function gi{x,Q'^). 
One of them is the Adler-Bardeen (AB) scheme which was introduced by requiring 
that the change from the MS scheme to this scheme be independent of x, so that 
the large and small x behavior of the gluon coefficient function is unchanged. In our 
case, we have in moment space 

Awin, a = AB) = , Az(n, a = AB) = Az(n, a = AB) = A^.- . (3.30) 

n 

(iii) [The off-shell (OS) scheme] In this scheme [^ we renormalize operators 
while keeping the incoming particle off-shell, ^ 0, so that at renormalization 
(factorization) point = — p^, the finite terms vanish. This is exactly the same as 
"the momentum subtraction scheme" which was used some time ago to calculate, 
for instance, the polarized quark and gluon coefficient functions [^, 30|. The CI- 



relations in Eqs.( 3.2^ ) and ( p.28| ) also hold in the OS scheme [51|, since the axial 



anomaly appears as a finite term in the calculation of the triangle graph for 
between external gluons (photons) and the finite term is thrown away in this scheme. 
The transformation from the MS scheme to the OS scheme is made by choosing 

Aw(n,a = OS) = CFl\Siin)f + 382(71) - S,in)(- - 

IL J [n + l) 
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Az{n,a = OS) 



7 2 3 

-o + 7 

2 n n + 1 

Az{n,a = OS) 
iV/(-4^5i(n) + ^ 



1 

— + 



+ 1) 



n 



(n + l)2 
4 

(n + 1) 



(3.3i; 



It is noted that in the OS scheme we have AA'^^ qs — ^^ns, os = o'^^y 
n = 1 but also for all n. 

(iv) [The Altarelli-Ross (AR) scheme] Using massive quark as a regulator for 



collinear divergence, Altarelli and Ross pT], |53] derived the same one-loop gluon 



coefficient function ACq as in the case of CI scheme. In order to obtain the one- 
loop quark coefficient function in this scheme, however, we need to do an extra 
subtraction so that the conservation of the nonsinglet axial currents is secured ||53| . 
The transformation rule is 



Awsin, a = AR) = Cpl 2 Si{n 



2S,in) 



1 1 

2 -| I 

n n + 1 {n + ly 



Si{n)[- 
2 

+ 



+ 2 



n 



Az{n, a = AR) = Az{n, a = AR) = 2A^^ 



3.4 The DIS^ scheme 



n(n + 1) 



(3.32) 
(3.33) 



An interesting factorization scheme, which is called DIS-y, was introduced some 
time ago into the NLO analysis of the unpolarized real photon structure function 
F2{x, Q^). Gliick, Reya and Vogt |39] observed that, in the MS scheme, the ln(l — x) 



term in the photonic coefficient function C^ix) for F2 , which becomes negative 
and divergent for x — > 1, drives the 'pointlike' part of F2 to large negative values 
as X — > 1, leading to a strong difference between the LO and the NLO results 
for pointlike largc-x region. They introduced the DIS^ scheme in which 

the photonic coefficient function C], i.e., the direct-photon contribution to F^, is 
absorbed into the photonic quark distributions. It is noted that, for the real photon 
target, the structure function F2 is decomposed into a 'pointlike' and a 'hadronic' 
part, the former being perturbatively calculable but not the latter. And beyond the 
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LO both the 'pointhke' and the 'hadronic' parts depend on the factorization scheme 
employed. A similar situation occurs in the polarized case, and the DIS^ scheme 
was applied to the NLO analysis for the spin-dependent structure function gi{x, Q"^) 
of the real photon target by Stratmann and Vogelsang |jl8 . 
In the polarized version of DIS^ scheme we take 

^ws{n, DIS^) = ^WNs{n, DIS^) = ^z{n, DIS^) = , (3.34) 

Az(n,DIS,) = ^Ai?;^ 



I n[n + 1) n n + 1 



so that 



/ 

= . (3.36) 

Note that the relation a la Eqs.(|3^) and (|3.27|) in the MS and Cl-like factorization 
schemes does not hold anymore in this scheme, i.e., 

Dis. ^ I^Ai^S, Dis. (= ^^Bl ^) . (3.37) 

For 77, = 1, we have 

Ai(n = 1,DIS^) = , (3.38) 

and thus, together with Eq. (|3.34| ), we observe that as far as the first moments are 
concerned, DIS^ scheme gives the same results with MS. In other words, in the DIS^ 
scheme, both the QCD and QED axial anomaly effects are retained in the quark 
distributions. 

With these preparations, we now examine the factoraization scheme dependence 
of the polarized parton distributions in the virtual photon. The two-loop anomalous 
dimensions of the spin-dependent operators and one-loop photon matrix elements 
of the hadronic operators in the MS scheme are already known. Corresponding 
quantities in a particular scheme are obtained through the transformation rules 
given in Eq.( p.l9| ). Inserting these quantities into the formulas given in Appendix 
A, we get the NLO predictions for the moments of polarized parton distributions in 
a particular factorization scheme. 



13 



3.5 Gluon distribution in the virtual photon 

Let us start with the gluon distribution. We find that all the factorization schemes 
which we consider in this paper predict the same behavior for the gluon distribution 
up to NLO: 

AG^(n, Q\ 1^ = AG^(n, Q\ P') |^ , (3.39) 

where a means factorization schemes of CI, AB, OS, AR and DIS^. This can be 
seen from the direct calculation or from the notion that, up to NLO, AC^la satisfies 
the same evolution equation as AC^Ij^ with the same initial condition at = P^, 
namely, AG^n, P^ P^)\a = AG''{n, P\ p2)|_ = q. 

If we consider a more general factorization scheme in which the hadronic part of 
Z^^in, Q"^) in Eq.( |3.8| ) is replaced with a new one as follows. 



1 ; V -f^Au 1 - ^Av 



(3.40) 



then, in this new factorization scheme, the predicted gluon distribution is not the 
same with AG"'{n, Q^, -P^)|ms NLO. However, the first moment is found to be still 
the same. In other words, the first moment of the gluon distribution in the virtual 
photon, AG'*'(n = 1,Q'^,P'^), is factorization-scheme independent up to NLO. This 
is due to the fact that the new terms, which appear by the inclusion of Am and Av, 



will be proportional to AK^'^ and that AK^'^~^ = 0. Also inclusion of Am and Av 
terms in Z^^ does not modify the photon structure function gKx, Q^, P^) itself up 
to NLO, since the gluon coefficient function starts in the order as- Moreover, the 
quark distributions in the virtual photon do not change by the inclusion of Am and 
Am terms. 



4 The n = 1 moments of parton distributions 

The first moments of polarized parton distributions in the virtual photon are partic- 
ularly interesting since they have relevance to the QCD and QED axial anomalies. 
The explicit expressions for the moments of Ag^, AG''', and Ag^^ up to NLO are 
given in Appendix A. We take the n — 1 limit in these expressions. Useful n = 1 
moments of one- and two-loop anomalous dimensions, photon matrix elements, and 
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coefficient functions both in the MS and Cl-hke schemes are enumerated in Ap- 
pendix B. As far as the ffist moments are concerned, DIS-y scheme gives the same 
results with MS. Note that we have 



0, Ar^ = -2/3o, Xl=s\^A^Tsn-0, (4.1) 

where A"^^ are eigenvalues of the one-loop hadronic anomalous dimension matrix 
A7°;""^. The zero eigenvalues A!p^ = — correspond to the conservation of 



the axial- vector current at one-loop order. 



4.1 The n = 1 moment of gluon distribution 

The expressions for the moments of gluon distribution are given in Appendix A. 2. 
In these expressions the factors 

11 1 



A!l ' 2/3o + XI ' 2Po + X1-X^ 



(4.2) 



may develop singularities at n = 1 and so we need a little care when we deal with 
them. Taking the limit of n going to 1, we find 







'G 



finite 



finite , ^ , B^'' finite , (4.3) 



Aa" 72{e^)NfCF ■ 

The terms proportional to L^", -B^", and Aq"' all vanish in the n — 1 limit, since 
they are multiplied by the following vanishing factors: 



1 - 



A!!/2/3o+l 



1 - 



Only exception is the term proportional to Aq^. We find for n 



2.1X^/2130 



Thus we obtain 



AG'^(n = l,g^p' 



12a 



a,m-c^s{P') 



(4.4) 



(4.5) 



(4.6) 



for the ffist moment of the gluon distribution in the virtual photon. It should be 
emphasized that the result is factorization-scheme independent. 
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4.2 The n = 1 moment of quark distributions 

The expressions for the moments of quark distributions are given in Appendix A.l 
and A. 3. In all the factorization schemes under study, i.e., MS, DIS^, and Cl-like 
schemes, we find for n — > 1, 



T +n 



Bp 



0, 

) 

, B 



finite , Ag 



S 



, 
" ^0 
finite 







A 



NS 



finite 




(4.7) 



The terms proportional to A^^" and B^^'' are multiplied by the vanishing factors in 
Eq.(|4.4|), and the A^g term multiplied by 



1 



A^c/2/3o 



and, therefore, the L^", L^'^, L%s, A^, A^", A^^, 5^", and B^g terms in 



Eqs.( |A.lD and ( [A.17| ) all vanish in the n = 1 limit. Then the first moments of quark 
distributions are given by 



a 



inpo 4:71 ^ 



(4.9) 
(4.10) 



We now see that scheme dependence for the first moments of quark distributions is 
coming from the photon matrix elements AA^ and AA^ 



1 n 



In the case of Cl-like factorization schemes, a = CI, AB, OS, AR, we have 
Aw{n = 1, a) = 0, 

Az{n = l,a) = Az{n = l,a) = Nf for a = CI, AB, OS, AR (4.11) 
We find from Eqs.( p.22D and ( p.25| ) that these schemes give 



(4.12) 
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This leads to an interesting result: The first moment of spin-dependent quark dis- 
tributions in the virtual photon vanish in NLO for a = CI, AB, OS, AR. 

AqUn = 1, Q^ P')\a = Aqlsin = 1, Q^ P^) 1 = (4.13) 

The vanishing first moments imply that the axial anomaly effects do not reside in the 
quark distributions. In these Cl-like schemes, the QCD and QED axial anomalies 
are transfered to the gluon and photon coefficient functions, respectively, and their 
first moments do not vanish. Indeed we obtain from Eqs.( p.l|) and 



ACS: T' 



AC 



7, n=l 
7, a 



71 



where we have used the fact [^ , 
1 



Nf 
f ^ 

1,12 



< > 



AC7^' 



MS 



MS 



TT 



A ^7. n=l 



(4.14) 

for a = CI,AB,OS,AR (4.15) 



MS 
2\ 



1 



TT 



0{al 



+ O{aai) 



On the other hand, in the MS (and also in DIS^) we obtain from Eq.i 

1 71=1 

V, MS 



and thus Aq^' " ^ 



A /1"=1 

NS, MS 

MS and AqJis=^ 



-12((e^) - {ey)Nf 



(4.16) 
(4.17) 

) 

(4.18) 
(4.19) 



MS 



are non-zero constant. Actually we can go 



one step further to the order of as QCD corrections. This is due to the fact that, in 
the MS scheme, the parton distribution Aq^in = 1)|ms = i^Qsy AC^, Aq1ig)\jj^ 
satisfies a homogeneous differential equation without inhomogeneous LO and NLO 
AK terms. Indeed we find 



AqUn=l,Q',P')\j^ 



Aq 



NS 



n 



TT 

a 



2 a,(P^) - a^iQ' 



-Nf U.20) 



^,Q',P')\ms = [-^3{{e')-{er)Nf]{l + 0{al)} , (4.21) 

the derivation of which is shown in Appendix C. In the MS scheme, the axial 
anomaly effects are retained in the quark distributions. The factors — ^ 3{e'^)Nf 
and — ^ 3(^(e^) — (e^)^^A^/ are related to the QED axial anomaly and a term 

is coming from the QCD axial anomaly. 



2 as{P')-as{Q^) 
/3o 



Nf in AgJ "='| 
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4.3 The n = 1 moment of gi{x, Q^, P^) 

The polarized structure function gi{x,Q'^,P'^) of the virtual photon satisfies the 



following sum rule||2y, |23 



f'dxgJ{x,Q\P' 
Jo 



■-{e'}N, ( 1 

TT V 



TT 



6a 



2 as{P') ~ as{Q') 



TT 



+ (4.22) 



This sum rule is of course the factorization-scheme independent. Now we examine 
how the scheme-dependent parton distributions contribute to this sum rule. In the 
Cl-like schemes (a = CI, AB, OS, AR), the first moment of the quark distributions 
vanish in NLO, and thus the contribution to the sum rule comes from the gluon and 
photon distributions. Equations ( [4.6|) and ( [4.15| ) show that 

l,Q',P')\a + AC2'r' 



'7, a 



(4.23) 



leads to the result ([4.22|) . On the other hand, in the MS scheme (and also in DIS. 



7/' 



the one-loop gluon and photon coefficient functions vanish, Ai?"; i-r^ 



Ai?"=J— = 

7, MS 



and, therefore, the sum rule is derived from the quark contributions. Indeed we find 
from Eqs.(|4l^), ( ^:20l - F2T| ) 

11, n=l rSl o2m I \n^,n=\ ^ „7 (, ^2 D2^ 



n 



IMS 



(4.24) 



leads to the same result. 

It is interesting to note that the sum rule ( |4.22| ) is the consequence of the QCD 
and QED axial anomalies and that in the Cl-like schemes the anomaly effect resides 
in the gluon contribution while, in MS, in the quark contributions. Furthermore, 
the first term of the sum rule ( [4.22| ) is coming from the QED axial anomaly and the 
second is from the QCD axial anomaly^. 



5 Behaviors of parton distributions near x = 1 

The behaviors of parton distributions near x = 1 are governed by the large-n limit of 
those moments. In the leading order, parton distributions are factorization-scheme 
^This notion was first pointed out by the authors of Ref.pO|| 
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independent. For large n, Ag2(n, Q^, -P^)|lo and AqJfg{n,Q'^, P'^)\i,o behave as 
l/(n In n), while AG^ {n, Q"^ , P'^)\i^o oc l/(n In n)^. Thus in x space, the parton 
distributions vanish for a; — > 1. In fact, we find 

-1 



AG\x,Q\P^)\^o 



"1 
a 



4 In (1 - x) 
—In X 



(5.1) 
(5.2) 



47ra,(g2)^-^^- '2 In" {I - x) ' 

The behaviors of AqJ^glx, Q^, P^) for x 1, both in LO and NLO, are always given 
by the corresponding expressions for Aqg{x, Q^, P^) with replacement of the charge 
factor (e^) with ((e^) - (eY). 

In the MS scheme, the moments of the NLO parton distributions are written in 
large n limit as 



AG\n,Q',P')Uo,m 
So we have near x — 1 

Agg(a;, Q^, -P^)Inlo, MS ~ 
AG^(x, (5^, P^)|nlo, MS ~ 



a , o, „ In n 
Att ^ n 



a 
An 



Nf{e^)6 -ln(l-a;) 



—In X 



(5.3) 
(5.4) 

(5.5) 
(5.6) 



It is remarkable that, in the MS scheme, quark parton distributions, Aq],{x)\^-^Q ^as 
and Aq]^g{x)\j^i^Q jj^ positively diverge as [— ln(l — x)] for x ^ 1. Recall that 
AG''{x,Q\P^)\ 

NLO is the same among the schemes which we consider in this paper. 
The NLO quark distributions in the CI, AB, AR and DIS^ schemes also diverge as 
a; — > 1, since their moments behave as In n/n in the large n-limit. We find for large 

X, 

a 



AqUx,Q',P')Uo, CI 
Ag2(a;, Q^P2) Inlo, ab 
Ag2(x, g^p2) Inlo, ar 
Ag^(a;,Q^P')|NLO, Dis, 



Nf{e')6 -ln(l-x) 



Ati 

a 

Ati 
a 

An 

^Nf{e')eHl-x) 



Nf{e^)Q -ln(l-a;)+2 



Nf{e^)lS -ln(l-x) 



(5.7) 
(5.8) 
(5.9) 
(5.10) 
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It is noted that Ag5(x, Q^, -P^)|nlo, dis^ negatively diverges as x ^ 1. This is due 
to the fact that the photonic coefficient function AC^{x), which in MS becomes 
negative and divergent for x — 1 , is absorbed into the quark distributions in the 
DIS^ scheme. 

On the other hand, the OS scheme gives quite different behaviors near x = 1 for 
the quark distributions. Since the typical two-loop anomalous dimensions in the OS 
scheme behave in the large n-limit as 



^INS, OS 

while in the MS scheme 

fl),n 



^7^! OS ^ n 



^^i-l;." — ^ Ao'^^'''^^^ — cx; In ri 

'nS, MS 'qq, MS 



n 



^^S, MS ^ 



In n 



n 



(5.11) 



(5.12) 



we find that the moment of Aq1{n, Q^, P^)|nlo in the OS scheme is expressed in the 
large rz-limit as 



AqUn,Q\P')\^^o, OS 



a 



69 



+ 7^/ 
4 ■ 



(5.13) 



In X space, Ag^(a;, Q^, P^)|nlo, os does not diverge for x — * 1 but approaches a 
constant value: 



Ag2(a;, Q^P') Inlo, os 



a 
47r 



69 



4 ■ 



(5.14) 



Therefore, as far as the large x-behaviors of quark distributions, and gluon and pho- 
ton coefficient functions (see Eqs.( |5.18| - 15.19| ) below) are concerned, the OS scheme 
is more appropriate than other schemes in the sense that they remain finite. Also 
the quark coefficient function in the OS scheme has a milder divergence for x ^ 1 
than those predicted in other schemes (see Eq.( ^.17|) ). 

Before ending this section, we now show that, as x — > 1, the polarized virtual 
photon structure function gj {x , Q'^ , P"^) approaches a constant value 



K 



a 

4:71 



51 3^ • 

h -Nf 

8 4^ 



(5.15) 



in NLO. The result is, of course, factorization-scheme independent. It is interesting 
to note that the constant value k coincides exactly with the one given in Eq.(4.39) 
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of Ref.[^], which was derived as the large n hmit of the moment of the NLO term 
62(3^) for the unpolarized structure function In the leading order, Eq.( [5.1| ) 

tells us that 



gUx, g^ Ilo = {e')Aql{x, Q^ P^) |lo + Aqlsix, Q^ P^) |lo 

a Atc , 4, 9 — 1 , , 

and thus gi{x, Q^, P^)|lo vanishes as x — > 1. 

In order to analyze the large x-bahavior of the next-leading order gj^x, Q^, P^) |nlc 
we need information on the coefficient functions. Note that AB^g\a = ABg\a. They 
behave, as a; — >• 1, 



AB.ix 



2Cf 
-2Cf 
3Cf 



2 lii(l-x) 
1-x 

2 In(l-x) 

1-x 



for a = MS, CI, AB, DIS^ , 
for a = AR , 
for a = OS , 



(5.17) 



ABGix)\a 



2Nf ln(l - x) for a = MS, CI, AB, AR, DIS^ , 



(5.18) 



-4Nf 



for a = OS 



AC{xyja 



^(e^)12Ar/ In(l-x) for a = MS, CI, AB, AR 







for a = OS 



for a = DIS^ . 



(5.19) 



The coefficient functions ABs{x) and ABg{x) are the same that appear in the po- 
larized nucleon structure function gilxjQ"^). The ABs{x) in all schemes considered 
here diverg 1, but OS scheme gives a milder divergence for ABs{x) than 

other schemes. Also note that ABg{x)\os remains finite as x 1, but ABg{x) in 
other schemes negatively diverge. 

I NLO 



Let us write gji^x, Q^, P^)|nlo in terms of partonic contributions as follows: 
gUx, Q\ P^)|nlo = W(x)|?,ro + 9l{x)\t^o + Aq(^) 



(5.20) 
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where 



i/i I-'-JInlo 



(e^) Ag2(x, Q^ P^) |nlo + Ag],^(x, Q\ P') |nlo 



LO 



47r 



An 



Aqls{x,Q\P')\LO 



-ABg{x)(^AG^{x,Q\P')\^o 



(5.21) 
(5.22) 



Then we find, for x 



i/i I-'-JInlo 



^(e4)12iV^ln(l 



x] 



8 ^ 4^ J 



for a 
for a 
for a 



MS, CI, AB, AR 
OS , 
DIS. 



(5.23) 



The NLO gluon contribution (^i (a;, Q^, P^)|nl^ vanishes faster than (In x)^ in any 
scheme under consideration. As for the NLO quark contribution, gKx, Q^, P^)|nlo^ 
in MS, CI, AB, AR schemes, diverges as [— ln(l — x)] for x ^ 1. However, Eq.( 5.19|) 
shows that the one- loop photon coefficient function AC^{x) in these schemes also 
diverges as [— ln(l — x)] with the opposite sign and the sum becomes finite. On the 
other hand, in the OS scheme, we observe from Eqs. (|5.23|) and ( p. 19 ) that both the 
quark contribution and photon coefficient function remain finite as x — »• 1, and it is 
easily seen that the sum 



5^1 (^)Inlo, OS + AC'(x)?^ OS 



(5.24) 



approaches the constant value k given in Eq. (|5.15| ). In the DIS^ scheme, the 
NLO quark contribution (7i (x)|nlo^ dis^ reaches the finite value k as x — 1, since 
Dis-y =0- fa.ct, as we see from Eq.( ^.21| ), gi{x)\'^{^Q is made up of two 
parts, the one from Ag^(x, Q^, P^)|nlo and the other from ABs{x)^Aq'g{x, Q^, P^)|lo5 
plus their non-singlet quark counterparts. In DIS^, both contributions diverge as 
X —>■ 1, but with the opposite sign, and the sum remains finite. 

The constant value k in Eq.( p.l5D is negative unless Nf > 9. Consequently, 
it seems superficially that QCD with 8 fiavors or less predicts that the structure 
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function gj^x, Q^, P^) turns out to be negative for x very close to 1, since the leading 
term gi{x, Q^, P^)|lo vanishes as x ^ 1. But the fact is that x cannot reach exactly 
one. The constraint {p + q)"^ > gives 

X < Xmax = g2 _^ p2 ' (5.25) 

and we find 

gUx = x^..,Q',P')\Lo > ^iV/(e')7^/?o (5.26) 
and the sum gj^x = Xmax, Q^, -P^)|lo+nlo is indeed positive. 

6 Numerical analysis 

The parton distribution functions are recovered from the moments by the inverse 
Mellin transformation. In Fig. 2 we plot the factorization scheme dependence of the 
singlet quark distribution Aq^i^x, Q^, P^) beyond the LO in units of 
{3Nf{e^)a/7r)\n{QyP^). We have taken Nf = 3, = 30 GeV^ P^ = 1 GeV^ 
and the QCD scale parameter A = 0.2 GeV. All four Cl-like (i.e., CI, AB, OS and 
AR) curves cross the x-axis nearly at the same point, just below x = 0.5, while 
the MS curve crosses at above x = 0.5. This is understandable since we saw from 
Eqs.( |4.13| , [4.20|) that the first moment of Ag^ vanishes in the Cl-like schemes while 



it is negative in the MS scheme. The DIS^ curve crosses the x-axis below x = 0.5, 
though the first moment of Ag^l^is^ is negative, taking the same value with the one 
in the MS scheme. Comparing the DIS^ curve at large x with the MS one, we will 
see that rapid dropping of the DIS^ curve as x — > 1 drives the crossing point below 
X = 0.5. 

As X 1, we observe that the MS, CI, AB, and AR curves continue to increase. 
In fact we see that the MS and CI curves tend to merge, the AB curve comes above 
those two curves and the AR curve diverges more rapidly than the other three. On 
the other hand, the OS and DIS^ curves start to drop at large x. The OS curve 
continues to increase till near x = 1, and then starts to drop to reach a finite positive 
value. The DIS^ curve reaches maximum at x ~ 0.8 and drops to negative values. 
These behaviors are inferred from Eqs.( p.5| , [5.7| - [5.10| , p.l4p . 
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Concerning the non-singlet quark distribution Aq'lfg{x,Q'^,P'^), we find that 
when we take into account the charge factors, it falls on the singlet quark distribu- 
tion in almost all x region; namely two "normalized" distributions Aq^ = Aq^/ (e^) 
and Aqjfg = AqJ^^/ {{e^) — (e^)^) mostly overlap except at very small x region. The 
situation is the same in all factorization schemes we have studied in this paper. This 
is attributable to the fact that once the charge factors are taken into account, the 
evolution equations for both Aq^ and AqJ^g have the same inhomogeneous LO and 
NLO AK terms and the same initial conditions at = (see Eq. ( |3.20D ). 

In Fig. 3 we plot again the OS and DIS^ predictions for Ag^(x, Q^, P^) together 
with the LO result. The motivation of having introduced DIS^ scheme into the 
analysis of the unpolaorized (polarized) real photon structure function (gj) was 
to reduce the discrepancies at large-x region between the LO and the NLO results 
for the 'pointlike' part of F2 {gj)- When applied to the polarized virtual photon 
case, it is seen from Fig. 2 and 3 that DIS^ scheme gives a better behavior for 
Aqg{x, Q^, P^) at large x than MS in the sense that DIS^ curve is closer to the LO 
result. However, we observe that absorbing the photonic coefficient function AC^ 
into the quark distributions in the DIS^ scheme has too much effect on their large- a; 
behaviors: The DIS^ curve for Aq'g{x, Q^, P^) goes under the LO one at x ~ 0.6 
and the difference between the two grows as x — > 1. In fact the DIS^ curve drops to 
negative values near at x = 1. 

From the viewpoint of 'perturbative stabilities' we find that the OS curve shows 
more appropriate behavior than the others. We see from Fig. 3 that the differences 
between the OS and LO curves are very small for the range 0.05 < x < 0.7. And 
the OS curve comes above the LO for x > 0.7. 

Fig. 4 shows the Q^- dependence of Ag^(x, Q^, -P^) in the OS scheme in units 
of (3A^/(e2)a/7r)ln(QVp2). Three curves with = 30, 50 and 100 GeV^ almost 
overlap in whole x region except in the vicinity of x = 1. We see from Fig. 4 that, 
in the OS scheme, Aq^ beyond the LO behaves approximately as the one obtained 
from the box (tree) diagram calculation, 

Ag^(^°^)(x,Q^p2) = (2x - l)3iV,(e^)^ In |^ . (6.1) 
The gluon distribution AG'^(x, Q^, P^) beyond the LO is shown in Fig. 5 in units 
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of (3A/'/(e2)a/7r)ln(Q7^^)> with three different values. Recall that every scheme 
considered in this paper predicts the same behavior for the gluon distribution up to 
NLO. We do not see much difference in three curves with different Q^. This means 
the AG"' is approximately proportional to ln(Q^/P^) . But, compared with quark 
distributions, AC is very much small in absolute value except at the small x region. 

In Fig. 6 we plot the virtual photon structure function gi{x, Q^,P^) in the NLO 
for Nf = 3,Q^ = 30 GeV^ and = 1 GeV^ and the QCD scale parameter A = 0.2 
GeV. The vertical axis corresponds to 

gUx.Q^P')/-Nf <e'>ln%. (6.2) 
Also shown are the LO result, the Box (tree) diagram contribution, 

^?7(^°^)(X,Q^P2) = {2x-l)-Nj <e'>\n%, (6.3) 

TT P^ 

and the Box diagram contribution including non-leading (NL) correction with mass 
being ignored 



TT 



(2x- l)ln^ -2(2x- l)(lna; + 1) 



(6.4) 

In our previous paper , there was an error in the program for numerical evaluation 
of the NLO gi{x,Q'^,P^). The corrected graph (NLO curve) here is different from 
the corresponding one in Fig.2 of Ref . ||23|| . The new NLO curve appears lower than 
the previous one for x < 0.7 and rather enhanced above x = 0.7. We observe 
that the corrected NLO curve remains below the LO one, and that the NLO QCD 
corrections are significant at large x as well as at low x. 

For the case of the real photon target, P^ = 0, the structure function can be 
decomposed as 

gj{x, Q^) = gj{x, Q^) Ipcrt. + gU^^ Q'^) Inon-pert. ■ (6.5) 

The first term, the point-like piece, can be calculated in a perturbative method. 
Actually, it can be obtained by setting P^ = in the expressions of parton dis- 
tributions in Eg. (|2.2|) or (|2.3|). The second term can only be computed by some 
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non-perturbative methods. In Fig. 7, we plot the point-hke piece of the real photon 
gi{x, Q"^) in the NLO, together with the LO result and the Box (tree) diagram con- 
tribution. The NLO curve, which is calculated by the corrected computer program, 
is different from the previous one in Fig. 6 in ref. The new NLO curve appears 
lower than the previous one for x < 0.6 and enhanced above x = 0.6. Also it remains 
below the LO curve. The NLO result qualitatively consistent with the analysis by 



Stratmann and Vogelsang |]T8|. In the unpolarized case, the moment of F2 has a 
singularity at n = 2 which leads to the negative structure function at low x. Thus 
we need some regularization prescription to recover positive structure function as 
discussed in Refs. ||5^, ^ ^ . Note that we do not have such complication at n = 1 
for the polarized case. 

Finally, in our numerical analyis, we took = IGeV^, which may not be 
necessarily large enough for the non-perturbative effects to be dying away. For our 
normalized parton distributions, however, the larger values of would not give any 
sizable change in shape and magnitude. 



7 Conclusion 

In the present paper, we have studied in detail the spin-dependent parton distri- 
butions inside the virtual photon, which can be predicted entirely up to NLO in 
perturbative QCD. The virtual photon target provides a good testing ground for 
examining the factorization scheme dependence of the quark and gluon distribu- 
tions. We have investigated the polarized parton distributions in several different 
factorization schemes. We derived the explicit transformation rules from one scheme 
to another for the coefficient functions, the finite photon matrix elements and the 
two-loop anomalous dimensions or parton splitting functions. 

In particular, we studied the QCD and QED axial anomaly effects on the first 
moments of quark distributions to see the interplay between the axial anomalies 
and factorization schemes. We find that, in the Cl-like schemes, the first moments 
of polarized quark distributions, both fiavor singlet and non-singlet, vanish in NLO 
while the standard MS scheme gives the non-zero value. Also we find that the large 
x-behaviors of polarized quark distributons dramatically vary from one factorization 
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scheme to another. Indeed, for x ^ 1, the quark distributions positively diverge or 
negatively diverge or remain finite, depending on factorization schemes. The numeri- 
cal analyses performed for the parton distributions reassures the above observations. 
From the viewpoint of 'perturbative stabilities' the OS scheme gives more appropri- 
ate behaviors for the quark distributions than the others. The gluon distribution 
turns out to be the same up to NLO among the six factorization schemes examined. 
Furthermore, its first moment is found to be factorization-scheme independent up 
to NLO. 

The same analysis on the factorization scheme dependence of the unpolarized 
parton distributions of the virtual photon can be carried out and will be discussed 
elsewhere. 
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Appendix 



A NLO expressions for polarized parton distribu- 
tions in the virtual photon 



We give the explicit expressions of Ag^, AG"'', and Ag^^ up to NLO. They are writ- 



ten in terms of one-(two-) loop anomalous dimensions Aj^-^ {A'-fl^^'"') =iIj,G), 



^Ins AKi'"" {AKi^^'"') {I = ■ip,G,NS), and the one-loop photon matrix 

elements of hadronic operators, AAf. The expressions of one-loop and MS scheme- 
two-loop anomalous dimensions are found, for example, in Appendix of Ref. . 



A.l Singlet quark distribution 
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(A.5) 
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X Ai^J" {(A7J; - 2/?o - A!!)A7S'" + A7^';A7;^^'"} (A7J; 
+AXj- {(A7J; - 2/3o - A!^)A7gi'" + A7S';A7g^"} A7j'S 
+2/3o(2/3o + A!^ - A!^) {A<)''^(A7J; - XI) + a4^)'"A7JS} 
-2/3o(2/3o + A^ - A^)A^AA;(A7j; - A^) 
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(2/5o + A!^)(A!^ 
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Xl){2Po + X 
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tplp 
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+ {(A7j;-A!!)A7 
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(A.9) 



A. 2 Gluon distribution 
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+ {(A7J; - A!^)A7S''^ + A7S;A7;;^-} A7S 
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(A.16) 



A. 3 Non- singlet quark 
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B The first moments 
B.l One- loop order 



where 



Mtt = A7j;=^ = (B.l) 

A7jr' = 0, A7S=^ = -6C^ (B.2) 
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A7Sr = -yC^+|r, = -2/3o (B.3) 

Xl=^ = 0, A!!=^ = -2/3o (B.4) 

AK^^T' = Aii'J"=^ = (B.5) 



4 A^f 
Ca = 3, C^-^^ Tf^^ (B.6) 



with A^j; being the number of flavors. 
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B.2 MS scheme 



- 24c^r, (B.8) 

= 8C^T, + ^ ^Cl = -2/3, (B.ll) 

. ^(l),n=l _ . ^(l),n=l _ . (l),n=l _ „ ,g 

^^;:hs = -12(e^)iV; (B.13) 

A^gTms = (B-14) 

A^FsIms = -me') - {er)Nf (B.15) 

Ai^rljs = ASr5:Ms = -3^i^ (B-16) 

^B-G:m - ^AS;=is = (B.17) 

(B.18) 

B.3 Cl-like schemes (CI, AB, OS, AR) 

AtS^T' = (B.19) 

= (B-20) 

A7irr^ = A7it-=0 (B.21) 

AT^rr - A7S':S = 18C^ - ^C^C^ + ^C^T, (B.22) 
AtS^::' = 32C^T^ + ^CaT^ - y Cl = -2A + 127V^C^ (B.23) 
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2 ^' " 
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(B.27) 
(B.28) 



C Derivation of Eqs.(S3g) and ( ^7211 ) 



We observe that, in the MS scheme, we have AK 



0,n=l 



AK" = (AK^, Ai^S, Ai^]^^). (Note AK^^%-[\^^ ^ 0. See Eq.(ill). ) Then, up 
to NLO, the parton distributions Aq^{n = 1)|ms ~ i'^ls^ AG^', Aq]^g)\j^ satisfy 
a homogeneous differential equation instead of an inhomogenious one: 



dAq^{n= 1,Q^P2) 



MS 



c/lng2 



Aq^{n = l,Q\P')\jj^AP{n=l,Q 



IMS 



(C.l) 



where the 3x3 sphtting function matrix AP is the hadronic part of AP given in 
Eq.(|3:T0|). Expanding AP{n = 1, Q^)|ms as 



APin=l,Q^)\^=^APi% + 



«'^W')l\p(i) |_ , 

^-'^n=llMS 
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and introducing t instead of as the evolution variable 



we find that Eq. (|C.l|) is rewritten as 

Aql^,{t)\^lAPi% + 
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/3o as{Q^y 

(0) , (^s{t) 



(C.2) 



(C.3) 



dt 
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A p(i) I Aap(°) 1 + Oio'^ 

^-r„=llMS OA? + '^'^"s 



(C.4) 
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We look for the solution in the following form: 



,7(1), 



with the initial condition (see Eq.(p.5| 
Agi°!(0) = 



(C.5) 



(C.6) 



a 

3a 



-^iV,((e^), 0, (e^)-(ey) (C.7) 



In the LO, we easily find that Aqr^^|(t) = due to the initial condition (|C.6|) . 
The evolution equation in the NLO is written as 



d Ag^ii (t) I MS _ A .7(1) /a p(0) ^ (t) 
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and we obtain for the solution 
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(C.8) 



(C.9) 



(C.IO) 



(C.ll) 



and using the information on the first moments of anomalous dimensions which are 
listed in Appendices B.l and B.2, we find that M turns out to be a triangular matrix 
in the following form: 



with 
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1 r as(O) _ asjt) 
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(l),n=£ 
Af5, MS 



TT TT 

The matrix exp(M) is, therefore, written in the form 

f e"" B 
exp(M) = e^^ 
V e"^ 

and thus we obtain from Eqs.( p.7| ) and ( |C.9| ), 
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where in the last hne we use the fact A^y^^^l^^ = 0. 
Incidentally, under the following approximtion. 



(C.19) 



b ^ -Cpt 



a + c ~ c 



(C.20) 



B is evaluated as 
B 



3Cf 



/3o 



(C.21) 



This leads to the expression for the first moment of gluon distribution AG^{n 
l,g2,p2^|MS given in Eq.il). 
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Figure Captions 

Fig. 1 Deep inelastic scattering on a polarized virtual photon in polarized e'^e" 
collision, e^e" — > e'^e~+ hadrons (quarks and gluons). The arrows indicate 
the polarizations of the , e" and virtual photons. The mass squared of the 
'probe' ('target') photon is -Q'^(-P^) (A^ < ^ Q^). 

Fig. 2 Factorization scheme dependence of the polarized singlet quark distribution 
Aql{x,Q'^,P^) up to NLO in units of {3Nf{e'^)a/7r)\n{Q^/P'^) with Nf = 3, 
g2 = 30 GeV^ P2 = 1 GeV^ and the QCD scale parameter A = 0.2 GeV, 
for MS (dash-dotted line), CI (solid line), AB (short-dashed line), OS (long- 
dashed line) , AR (dashed line) and DIS-y (dash-2dotted line) schemes. 

Fig. 3 The polarized singlet quark distribution Aqg{x, Q'^,P'^) up to NLO predicted 
by the OS and DIS^ schemes in units of {3Nf{e'^)a/n)\n{QyP'^) for Nf = 3, 
g2 = 30 GeV^ P2 ^ I Qgv^ and A = 0.2 GeV, together with the LO result. 

Fig. 4 The polarized singlet quark distribution Aqg^XjQ"^, P'^) up to NLO in the 
OS scheme in units of {3Nf{e^)a/7T)\n{Q'^/P'^) with three different Q'^ values, 
for AT^ = 3, P2 ^ 1 GeV^ and A = 0.2 GeV. 

Fig. 5 The polarized gluon distribution AG'^{x,Q'^,P'^) beyond the LO in units 
of {3Nf{e^)a/7v)\n{QyP^) with three different values, for Nf = 3, P^ ^ 
1 GeV^ and A = 0.2 GeV. 

Fig. 6 Polarized virtual photon structure function gi{x, Q"^ ^P"^) up to NLO in units 
of (3A^/a(e^)/7r) ln(QVp2) for = 30 GeV^, and P^ = 1 GeV^ and the QCD 
scale parameter A = 0.2 GeV with Nf = 3 (solid line). We also plot the LO 
result (long-dashed line), the Box (tree) diagram (2dash-dotted line) and the 
Box including non-leading contribution. Box (NL) (short-dashed line). 

Fig. 7 Point-like piece of the real photon structure function gi{x, Q"^) in NLO in 
units of {3Nfa{e^)/-K) ln(gVA2) for Q"^ = 30 GeV^ with A = 0.2 GeV, Nf = 3 
(solid line). Also plotted are the LO result (long-dashed line) and the Box 
(tree) diagram contribution (short-dashed line). 
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